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ABSTRACT 

We report the discovery of fast variability of 7-ray flares from blazar OJ 287. This 
blazar is known to be powered by binary system of supermassive black holes. The 
observed variability time scale T var < 3 — 10 hr is much shorter than the light crossing 
time of more massive (1.8 x 1O 1O M ) black hole and is comparable to the light crossing 
time of the less massive (1.3 x 10 8 M©) black hole. This indicates that 7-ray emission 
is produced by relativistic jet ejected by the black hole of smaller mass. Detection 
of 7-rays with energies in excess of 10 GeV during the fast variable flares constrains 
the Doppler factor of the jet to be larger than 4. Possibility of the study of orbital 
modulation of emission from relativistic jet makes OJ 287 a unique laboratory for the 
study of the mechanism(s) of formation of jets by black holes, in particular, of the 
response of the jet parameters to the changes of the parameters of the medium from 
which the black hole accretes and into which the jet expands. 
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1 INTRODUCTION 

Although large-scale jets ejected by Active Galactic Nu- 
clei (AGN) were discovered almost a century ago, origin of 
this phenomenon remains o bscure up to present days (see 
lHarris feKrawczvnskil (|2006l ) for a recent review). It might 
be that the jets are accelerated via magneto-centrifugal 
force along twisted mag netic field lines above accre tion disk 
around the black hole (jBlandford fc Pavnd Il982l ) . Other- 
wise, an outflow can be created via Blandford-Znajek mech- 
anism of electromagnetic power extraction from a rotating 
BH, similar to the mechanis m responsible for the gene ration 
of relativistic pulsar winds (|Blandford feZnaiekll 19771 ). 

Blazar OJ 287 (2 = 0.306 (|Stickel et all 1 19891 ) ) pro- 
vides a unique laboratory for the study of the mechanism 
of AGN activity, because this is on e of the few AGN known 
to host binary black hole system (|Lehto fc Valtonenlll996l : 
IValtonen et"aU 120091 ). In this system, a lighter black hole 
of the mass Mbhi — 1.3 x 10 s Mq orbits a heavier black 
hole of the ma ss Mbh2 — 1.8 x 10 10 M fl with a perio d 
Porb ~ 11.65 yr (|Valtonen et~aHl2009l : ISillanpaa et aO 1988h . 
Separation of the components of the system at periastron is 
just about 10 Schwarzschild radii of the heavier black hole, 
so that the orb ital motion is strongly affected by relativistic 
gravity effects (|Valtonen et al.|[2008f ) . 

OJ 287 is known to belong to the BL Lac sub-class of 
AGNi, which means that it emits a relativistic jet whose 
direction is aligned with the line of sight. It is not clear 



a priori, which of the black holes ejects the observed rel- 
ativistic jet. Following a naive argument, which does not 
take into account the relativistic beaming of the jet emis- 
sion, one would assume that the observed relativistic jet is 
the one ejected by the heavier black hole, simply because 
the bigger black hole accretes more matter and, therefore, 
could produce more powerful jet. Most of the existing stud- 
ies of multi-wavelen gth blazar activi t y of O J 287 adopt this 
assumption (see e.g. Valtonen et alj ( 20091 )). However, rela- 
tivistic jets in BL Lacs are known to Doppler factors i > 1. 
This results in boosting the apparent luminosity of the jets 
by a factor <5 4 , when the jets are viewed face-on. Thus, if the 
less powerful jet emitted by the smaller black hole is aligned 
with the line of sight, while the jet from the larger black 
hole is not, the jet from the smaller black hole might give 
dominant contribution to the source flux. 

In what follows we show that variability properties of 
7-ray emission from the source indicate that the relativis- 
tically beamed emission comes from the jet produced by 
the smaller black hole. Independently of the value of the 
Doppler factor of the jet S, the shortest observed variabil- 
ity time scale AT m j n imposes a constr aint on the size of 
the jet's "central en gine", R C e < cAT (jCelotti et al.lll998l : 
iNeronov et~al]|2008l ). In the case of OJ 287, R C e turns out 
to be much smaller than the Schwarzschild radius of the 
more massive black hole, but compatible with the size of 
the smaller mass black hole. 

Observation of 7-ray emission from the base of the jet of 
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Figure 2. Lightcurve of October-November 2009 flare of OJ 287 
in E > 0.1 GeV energy band, binned in time bins with S/N ratio 
equal to 3. Black curve shows model fit to the lightcurve, Eq. JTJ 
with parameters given in Table [T] 



the smaller black hole in the system makes OJ 287 a unique 
laboratory for the study of mechanisms of jet production. 
Regular orbital modulation of the physical parameters of 
the ambient medium around the 1.3 x 1O 8 M0 black hole 
provides a unique possibility to study the response of the 
jet to the changes of the properties of accretion flow and 
of external medium in which the jet propagates. In this re- 
spect, the OJ 287 system provides a scaled-up analog of 
Galactic 7-ray-loud binaries, in which orbital modulation 
of 7-ray emission enables a study of response of relativistic 
outflow from a compact object (a neutron star or a black 
hole) to the changes of the properties of external medium 
(stellar wind and radi ation field of companion star) (see e.g. 
IZdziarski et~aLl (j2010T l). 



2 FERMI OBSERVATIONS 

In order to study variability of the 7-ray signal during the 
flaring activity, we have processed publicly available data 
of the LAT instrument, using the Fermi Science Tools pro- 
vided by the Fermi Science Support Centre. The data were 
selected using gtselect tool. The lightcurves were produced 




Figure 3. Images of sky region around OJ 287. Panels a and 
c show 0.3-1 GeV and 1-300 GeV band images for the entire 
observation period. Panels b and tishow 0.3-1 GeV and 1-300 GeV 
band images during the flaring activity of OJ 287. Coordinate 
grid show ecliptic coordinates. Green circle shows the region from 
which the background lightcurve was extracted. Images in the 0.3- 
1 GeV band are smoothed with a gaussian of the width 1° . Images 
in the 1-300 GeV band are smoothed with the Gaussian with 0.5° 
width. 



with the help of gtbm and gtexposure tools as it is explained 
in the Fermi data Analysis Threads Q. 

Fig. [T]shows the longterm lightcurve of the source at the 
energies above 0.1 GeV over the August 2008 - September 
2010 period binned to achieve signal to noise ratio S/N=5 
per time bin. Photons were collected from the circle 2° in 
radius and centered on OJ 287 . The flaring activity period, 
which started on October 5 2009, could be readily identi- 
fied. The detailed lightcurve of the flaring period shown in 
Fig. [2] reveals several well separated flares. The brightest 
flare which h a ppene d on October 22, 2009 was reported by 
ICiprini et al.l (|2009). Followup observations of October 22 
flare by Swift/XHT have reveale d an increase of the X-ra y 
flux accompanying the 7-ray flare (|D'Ammando et al.ll2009l ). 

The direction toward OJ 287 is situated close to the 
plane of Ecliptic. The lightcurve of the source could be af- 
fected by the passage of the Sun through the field of view 
of the telescope. The periods of the Sun passage within the 
region of the radius 15° centered on OJ 287 (the sky region 
chosen for the analysis) are shown by the orange-shaded 
strips in Fig. [1] Possible effect of the Sun passage on the 
source lightcurve is most clearly seen in the lightcurve col- 
lected from a region of the radius 2° displaced by ~ 5° 
from the position of OJ 287 to the position RA=127.11 
DEC=19.89. No source is detected at this position, so the 
lightcurve collected from this region (shown by grey data 
points in Fig. [1} could be considered as a measure of the 
diffuse background level close to the position of OJ 287. One 
could notice a flare in the background lightcurve, associated 
with the Sun passage close to the background extraction re- 



1 http : ll± ermi . gsf c.nasa. gov/ssc/data/analysis/scitools/ 
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Table 1. Parameters of the model fit (see Eq. JTJ) to the 
lightcurve of the 7-ray flare of OJ 287 together with their 68% 
confidence ranges. 



gion. The effect of the Sun passage on the lightcurve of OJ 
287 is less pronounced (see Fig. [1]), because the source is 
situated higher above the Ecliptic plane, compared to the 
background region. 

Analysis of the images of 7-ray sky around OJ 287 show 
that, apart from OJ 287, the region of the radius 15° cen- 
tered on OJ 287 contains several other sources (most clearly 
visible in the 1-300 GeV band image in the panel c of Fig. 
[3|. However, these sources are rather weak, so that they are 
not detected on one month exposure time scale correspond- 
ing to the duration of the flaring activity of OJ 287(see the 
right panel of Fig. [3]). Absence of strong sources near OJ 287 
ensures that the observed variability of the signal from the 
source during the flaring period is not affected by possible 
variable emission from a nearby source for which the tail of 
the point spread function overlaps with that of OJ 287. 

The E > 0.1 GeV lightcurve of the flare consists of sev- 
eral well separated pulses with rather sharp rise and decay. 
To find the rise and decay times we have fitted the lightcurve 
with a phenomenological model of a sum of exponentially 
rising and decaying pulses 



F(t) = B + Y, 



A k exp ((i - tk) /Uk) , t<tk 
A k exp(-(t-tk)/t dk ) , t>t k 



(1) 



where B — const is the background level. The background 
level was found from the circle of the radius 2° displaced by 
5 degrees from the source position. Parameters of the model 
function (0, derived from the fitting, are given in Table [1] 
One can see from Fig. [2] and Table [1] that brightest 
flares are characterized by the rather sharp rises and de- 
cays, with the rise/decay times of several hours. At this 
time scales measurement of the rise/decay times of the flares 
is complicated by the fact that Fermi/LAT telescope ob- 
serves a given patch of the sky once in 3.2 hours (once 
per two rotation periods of ~ 96 min). This is clear from 
Fig. U in which lightcurve for the time interval of several 
hours around the brightest flare from the source is shown 
in more details. The upper panel of the figure shows the 
E ^ 0.1 GeV lightcurve of the source. The lower panel 
shows the energies of photons collected from the circle of 
the radius 2° centered on the source (black points) and 
photons collected from a background extraction circle dis- 
placed of the same radius displaced by 5 degrees from the 
source position. Photons from the source and background 
regions come only within periodic time intervals spaced by 
3.2 hr marked by vertical grey strips in the two panels 
of the Figure. It is clear that Fermi/LAT pointing pat- 
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Figure 4. Upper panel: lightcurve of brightest episode of 
October-November 2009 flare of OJ 287 in E > 0.1 GeV en- 
ergy band. Lower panel: energies and arrival times of 7-ray from 
the source (black points) and from the background region (grey 
points). Vertical grey strips show the periods when the source was 
in the field of view of LAT telescope. 



tern does not allow to constrain the rise/decay time of the 
flares to better than 3.2 hr. The apparently abrupt end of 
the flare (12 photons from the source are detected between 
t =MJD55126.0+3 hr and t =MJD55126.0+4 hr, while no 
photon is detected within the subsequent observation pe- 
riod MJD55126.0+6.2 hr< t <MJD55126.0+7.2 hr) limits 
the decay time of the flare to be less than 3.2 hr. Assuming 
that the source flux did not change in the time interval fol- 
lowing the peak of the flare, one could estimate the chance 
probability of detecting zero photons in this time bin to be 
7 x 10 -4 , taking into account the exposures in the two ad- 
jacent time intervals with maximal and zero count rates, 
3.3 x 10 6 cm 2 s and 2 x 10 6 cm 2 s, respectively. 

To summarize, both fitting of the lightcurve profile with 
a phenomenological exponential rise / exponential decay 
model and and direct photon counting in individual 1 hr 
long Fermi exposures indicate that the flux from the source 
is significantly variable on time scales shorter or compara- 
ble to 3-10 hr. This fact has important implications for the 
physical model of the origin of 7-ray emission from the OJ 
287 system. 



3 ORIGIN OF RELATIVISTIC JET IN OJ 287 

Constraint on the variability time scale of the 7-ray flares 
T var = min(t r k,tdk) ^ 3.2 hr, derived above, enables to iden- 
tify the 7-ray emission site within the binary black hole sys- 
tem of OJ 287. 

It is commonly accepted that 7-ray emitting jets are 
generated by the AGN "central engines", the supermassive 
black holes, on the distance scales of the order of gravita- 
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Figure 5. Comparison of the upper bound on variability time 
scale of OJ 287 flares (grey shaded region) with the light crossing 
times of the two supcrmassive black holes in the system. 
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where A/bh is the black hole mass. Minimal variability time 
scale of electromagnetic emission originating from the AGN 
central engine is expected to be not shorter than the light- 
crossing time of the supermassive black hole, 



Tic = (1 + z)2Rbu/c = 2(1 + z) ( R g + 



7? 2 -a 2 )/c(3) 



0.5 [Mbh/1.3 x 10 s M Q ] hr, 
0.9 "Mbh/1.3 x 10 s M Q ] hr, 



a = R g 
a = 



where _Rbh is the size of the black hole horizon and < 
a < R g is black hole rotation moment per unit mass. Vari- 
ability of X-ray emission at time scale T ~ T\ c is observed 
in X-ray emission from Galactic sources powered by black 
holes with masses Mbh ~ WMq |Remillard fc McClintockl 
l2006t ). Variability at the time scale T var ~ Ti c is ob- 
served also in 7-ray emission from blazars, a speci al type 
of AGN with jets aligned al o ng the line of sight ( Gaidos 
19961; Aharonian et alj I2OO6I; iNeronov fc Aharonianl 12007 



Albert et al. I 120071 ; lAharonian et al. I 120071 ; INeronov et al I 
200ct l. 

Fig. [S] shows a comparison of the upper limit T var ^ 
3.2 hr on the variability time scale of the flares from OJ 
287 with the light crossing time of the two supermassive 
black holes in the OJ 287 system. One can see that the 
upper bound is much lower than the light crossing time of 
the more massive black hole with Mbh — 1.8 x 10 Mq, 
independently of the black hole rotation moment a. At the 
same time, the limit on the variability time scale is larger 
than the light crossing time of the smaller black hole. This 
implies that the observed 7-ray emission is generated by the 
jet emitted by the lower mass companion black hole in the 
system. 

Note, that the above result does not depend on the as- 
sumptions about relativistic motion of the 7-ray emission 
region. Indeed, the minimal possible time scale AT m ;„ of 



emission from relativistic jet moving with bulk Lorentz fac- 
tor r does not depend on F and is instead determined by 
the siz e of the non-moving "central engine" w hich ejected 
the jet jCelotti et al.lll998l ; lNeronov et alJ|200Sl lah. In the jet 
comoving frame, the minimal variability time scale AT' min ~ 
rAT mm is the light-crossing time of the smallest size 7-ray 
emitting "blobs", A.T' mirl ~ R' min /c. 7-ray emitting blobs 
which were ejected by the jet's "central engine" of the size 
Rce could not have size smaller than R' m i n -RceF in the 
comoving frame. This implies that AT m i n ~ AT^^/F ~ 
R'min/cT ^ Rce/c, independently of T. 

At the first sight, such conclusion looks counter- 
intuitive. Indeed, since the two black holes accrete matter 
from the same reservoir, the accretion rate onto more mas- 
sive black hole should be much higher. In general, higher 
accretion rate should lead to production of a more powerful 
jet. 

However, the above argument does not take into ac- 
count possible relativistic beaming effects, which are nor- 
mally very significant in BL Lacs. 7-ray flux from relativis- 
tically beamed jet moving with bulk Lorentz factor V at an 
angle 9 with respect to the line of sight is boosted by a factor 
S 4 , where 5 = (T(l—/3 cos is the Doppler factor and /3 is 
the bulk velocity of the jet. A constraint on the bulk Lorentz 
factor of the jet produced by the Mbh — 1.3 x lO 8 Af0 black 
hole could be found under assu mption that the o bserved X- 
ray emission from the system (|Seta et al.l I2009T 1 originates 
from the same jet region as the GeV emission. 

The highest energy of 7-ray photons coming from the 
source is E li7n ax ^ 30 GeV (see Fig. 7-rays of such en- 
ergies can produce e + e~ pairs in interactions with X-ray 
photons of the energies 
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(4) 



where we have assumed that typical collision angles for pho- 
tons emitted from the jet are a ~ Apparent (relativis- 
tically beamed) luminosity of OJ 287 in the soft X-ray band 
is Lapp — 3 x 10 44 erg/s, whi ch corresponds to the observed 
flux F x ~ 10" 12 erg/cm 2 s (|Seta et all 120091 ). Calculating 
the optical depth of 7-ray emission region of comoving size 
R' ~ 5cT V ar/(l + z) and luminosity L' ~ 5~ 4 (1 + z) 2 Lx one 
finds 



0.5 



10 44 erg/s 



3.2 hr 



(5) 



High energy 7-rays could escape from the source if r 77 < 1. 
This condition imposes a restriction on the Doppler factor 
S > 4. 

Thus, the observed 7-ray flux from the jet is Doppler 
boosted by at least a factor 8 > 3 x 10 . It is interesting to 
note that even if emission from the higher mass black hole is 
not relativistically beamed toward observer on the Earth, it 
might be noticed in the spectrum of the source. Indeed, as- 
suming a simple Eddington-like scaling of the accretion rate 
and jet luminosity with the black hole mass, Li ~ Mbhi, one 
finds that the relativistically beamed luminosity of the jet 
from the lighter black hole L ~ 5 4 Li ^ 10 2 Li could, in fact, 
be comparable to the overall luminosity of the heavier black 
hole, L 2 ~ (Mbh2/M B hi)£i ~ 1.4 x 10 2 Li. Fast variability 
of emission could therefore, serve as a tool for identification 
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of the contribution of emission from the lighter black hole 
in the overall source spectrum. 

A common feature of all models of jet production by 
black holes is that matter ejection into the jet is associ- 
ated with rotation of matter aro und the black hole and /or 
with rotation of the blac k hole (|Blandford &Znaiek|[l977l ; 
iBlandford fc Pavnel fl982l ) . This implies that characteristic 
time scale at which the properties of the jet could change 
is given by the period of rotation of the black hole itself or 
of the accretion flow onto the black hole. Period of rotation 
around a ci rcular orbit at a dist ance r from the black hole 
is given by l|Bardeen et al.lll972h 



r 3/2 ±aR 1/2 
P(r) = 2*(1 + Z ) r ±a f" 

CRg' 



(6) 



The + (— ) sign corresponds to the prograde (retrograde) 
orbit. Stable circular orbits exist only down to certain dis- 
tance r ms from the BH. The period of rotation along the last 
prograde stable orbit at the distance r ma is 



f(r- ms ) 



Mbh 



22 



1.3 x 10°M Q 
Mbh 



hr. 



a = R a 



1.3 x lO^M© 



(7) 



hr, a = 



Upper bound on the variability time scale T var ^ 3.2 hr 
is much shorter than period of rotation around the non- 
rotating black hole and is comparable or smaller than the pe- 
riod of rotation around maximally rotating black hole with 
a = R g . This means that relativistic ejections into the jet, 
responsible for the observed flares, are produced by the mat- 
ter moving in the direct vicinity of the black hole horizon, 
well inside the R = 6R g radius of the last stable orbit around 
non-rotating black hole. 



4 CONCLUSIONS 

To summarize, we find that observations of OJ 287 in the 
E > 0.1 geV energy band constrain the minimal timescale 
of flux variations of the source to be shorter than 3.2 hr. The 
upper limit on the minimal variability timescale imposes a 
restriction on the size of the jet formation region. We find 
that the size of the jet formation region in the OJ 287 sys- 
tem is much smaller than the size of the horizon of the more 
massive black hole in the binary black hole system power- 
ing the source. This means that the observed 7-ray emission 
is produced by the jet ejected from the smaller mass black 
hole. Higher apparent luminosity of the smaller mass com- 
panion is explained by the effect of relativistic beaming of 
7-ray emission. Combining X-ray and 7-ray data, we find a 
restriction on the Doppler factor of the 7-ray emitting part 
of the jet, 5 > 4. The observed variability time scale indi- 
cates that relativistic jet from the smaller mass black hole 
is formed close to the black hole horizon, well inside the last 
stable orbit around non-rotating supermassive black hole. 

7-ray data provide a new insight in the physical model 
of the unique binary supermassive black hole system in OJ 
287. 7-ray flaring activity is produced in connection with the 
passage of the smaller mass black hole through the accretion 
disk around the larger mass black hole during the periods 
of close approach of the two black holes in the periastron 
of the binary orbit. Interaction of the smaller mass black 



hole with the larger mass black hole accretion disk leads to 
the transient episodes of ejection into relativistic jet from 
the smaller mass black hole. It appears that the transient 
jet from the small mass black hole happens to be aligned 
along the line of sight, the fact responsible for the BL Lac 
type appearance of the source. It is not clear a priori if the 
jet from the smaller mass black hole forms only during the 
periastron passage or it exists throughout the binary orbit. 
If the jet is powered by the transient accretion onto the black 
hole, it would be natural to expect that the the small mass 
black hole jet (and the associated 7-ray emission) should 
disappear soon after the periastron passage on characteristic 
time scale of accretion. Otherwise, if the jet is powered by 
the rotation energy stored in the small mass black hole, it 
is natural to expect that the jet and the 7-ray from the jet 
would be persistent throughout the binary orbit. Systematic 
monitoring of the source evolution in 7-rays on the orbital 
(11.7 years) time scale, which is now possible with Fermi, 
might clarify this question. It is interesting to note that if 
the jet from the smaller mass black hole is directed along the 
black hole spin axis, its alignment with the line of sight might 
be destroyed by the pr ecession of the black hole spin axis (see 
IValtonen et al.l |2006f ) for detailed discussion of the orbital 
evolution of the system). This would mean that the BL Lac 
appearance of of the source might be time-dependent. An 
immediate consequence of the mis-alignment of the smaller 
black hole jet with the line of sight should be the loss of 
the strong Doppler boosting of the flux. In the absence of 
Doppler boosting, the emission from the smaller mass black 
hole might become sub-dominant compared to the emission 
from the larger mass black hole. Study of the details of the 
overall time evolution and of the short time scale variability 
properties of the source along the orbit and from periastron 
to periastron should clarify the transient/permanent nature 
of the BL Lac appearance of the source. 
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